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Abstract 

We have studied a Hubbard ring with alternating site potentials for half-filling in 
presence of a magnetic flux. Using a mean-field approach we have calculated the con- 
ductivity of such a ring at low and high temperatures. The interplay of the correlation, 
the polarizing field and the chemical modulation in site potentials tune the conductivity 
in an interesting fashion. In presence of the modulation in the site energy an appre- 
ciable variation in the conductance is observed with the change in flux. This scenario 
gets substantially modified in presence of the Hubbard correlation. Finite-size effects are 
also identified and they are found to be quickly disappearing with increasing system size. 
Sharp changes in the magnetoconductance is found to disappear at higher temperatures. 



I. Introduction 

The role of an externally applied magnetic field on the electrical transport properties of 
a metal has been extensively studied for a long time [1]. Recent studies of unusual magne- 
toresistive properties of low-density Si MOSFET's [2] and p-GaAs/AlAs heterostructures 
[3] have generated a great deal of interest from both theoretical and experimental points of 
view [4]. Quasi-two-dimensional organic systems e.g. /3"-(BEDT-TTF)2SF5CH2CF2S03 
also show interesting field-dependence of the electrical transport properties [5]. A similar 
type of field-dependent phenomenon is observed very recently in highly oriented pyrolitic 
graphite [6] . Although a clear theoretical picture is still to emerge in this area the recent 
studies clearly point out that the electronic correlation is going to play a very important 
role in the magnetoresistivity of some of these dilute two-dimensional (2D) systems (e.g. 
Si MOSFET) [4]. 

On the other hand persistent currents in mesoscopic metallic rings have been well 
studied for a long time from both theoretical [7] and experimental [8] perspectives. Such 
a ring, penetrated by a magnetic fiux, is found to show a typical Aharonov-Bohm type of 
oscillation [9] in the current. The presence of disorder is expected to enrich the problem 
[10, 11, 12] due to impurity scattering. More recently the effect of electron-electron 
interaction [13, 14] on mesoscopic conductance has generated a great deal of interest. 
Also a huge activity is going on with the manganite perovskites (e.g. Lai-^jSr^; MnOs or 
Lai-xCaxMnOa ) for their noticeable magnetoresistive properties [15] known as colossal 
magnetoresistance. 

In view of the scenario nientioned above we investigate the effect of the magnetic 
field on the conductivity of a mesoscopic ring which otherwise shows strong effects of 
electronic correlation in its ground state. A Hubbard chain with a periodic variation of 
the site potential is found to undergo a spin density wave (SDW)/ charge density wave 
(CDW) transition at half-filling [16, 17], and the conductivity of such a chain of finite 
length captures interesting features of the competition between the SDW and the CDW 
[17] instabilities. The effect of an externally applied magnetic flux on the conductivity 
in such a ring is, therefore, expected to reveal some significant aspects of mesoscopic 
systems. 
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The present problem requires some suitable approximate scheme for solution. Here we 
have adopted a mean field approximation (namely, the generalized Hartree-Fock approx- 
imation, GHFA in short) to tackle the problem. It is to be mentioned here that a mean 
field approach was found to be very much useful for capturing the SDW/CDW transition 
in the present model in absence of the magnetic field [16, 17]. 

II. The Model and the GHFA calculations 

The chemically modulated Hubbard chain has alternate sites occupied by A and B 
types of atoms. We take a finite chain of A'' sites, where N is an even number, in the form 
of a ring which is threaded by an externally applied magnetic field H; the Hamiltonian 
looks like: 

i \ ia / i i 

where e, is the site energy at the i-th site; it takes on the value ea or es depending on 
whether it belongs to an odd (A-type) or an even {B type) numbered site respectively. 
The number operator 77.^0- = cj^Cia- (cr |) and rii = ni-f + n^^; t is the hopping integral 
between the nearest neighbour sites. The phase factor (e"^^) appears because of the flux 
{—N(f)i in units of basic flux quantum (pQ = hc/2TTe) threaded by the ring due to the 
applied magnetic field. U is the on-site Coulomb interaction. The fourth term represents 
the Zeeman interaction between the local moments (in units of Bohr magneton) and 
the applied magnetic field H {—ATT(f)i/N for unit lattice spacing). In order to study the 
conductivity of the present model we calculate the response to a small static electric field 
in terms of the Drude weight [18, 19]. This amounts to the application of an additional 
infinitesimally small flux 02 (in units of basic flux quantum 0o = hc/2Tie ), which modifies 
the phase factor in the hopping term as well the Zeeman term, and the Hamiltonian looks 
like : 

n^Y. ^i'^i + 1 U^'^'^f"^ E + h.c) + [/ ^ n^tna - /f(02)/2 ^ (n^^ - n,|) (2) 

i \ i(7 / i i 

On employing the GHFA the Hamiltonian is decoupled into up-spin and down-spin 
parts with the modified site energies for the up- and the down-spin electrons. The modified 
site energies are given by: 

e'if ^ei + U < riii > , e-| = -I- [/ < n^i > (3) 
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This amounts to a decoupling of the Hamiltonian in to down and up spin parts as foUows: 

= E(4 - ^/2)«2.T + El^A - H/2)n2^-n + E c\^c,+r^e^^ + h.c) 

i i i 

ni = ^{e's + H/2)n2ii + EC^a + H/2)n2i-n + (< E 4iCi+ne^'^ + h.c) (4) 

i i i 

where 4'=4'i + (f)2- Here e'g^ and e'^^ are determined according to equation (3). 

Upon diagonalization (in a self-consistent method) of the decoupled "independent 
electron" Hamiltonian (4) one can calculate the partition function 

where /? = ^/kT, k being Boltzmann's constant, and T, the absolute temperature, jj, is the 
chemical potential. This leads to direct calculation of several thermodynamic quantities. 
The Drude weight can be evaluated from 



N ^ 

n 



(5) 



Here -Erie's are the single particle energy levels obtained by diagonalizing the decoupled 
Hamiltonian (4) for the spin species a. 

The CDW order parameter C, the SDW order parameter S and the magnetization M 
are defined by: 

c = ^E(-irK+"u-i) (6) 

i 

5=^E(-imT-na) (7) 

i 

^ = ^E(«»T-^a) (8) 

i 

We now discuss the special case oiU — which is exactly solvable for the infinite system 
and also gives an insight into the finite size case, for both U — Q and C/ > 0. 

In this case {U—0) exact diagonalization of the Hamiltonian (4) yields two energy 
bands for both species of spins. The bands for up and down spin take the following 
shapes respectively: 

E{k)^ = {eA + eB)/2 - H/2 ± ^^(e^ - ea? + 8*2(1 + cos 2{k + 0)) 
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E{k)i = {eA + eB)/2 + H/2 ± ^{63 - 6^)^ + 8^2(1 + cos 2{k + 0)) (9) 

where k is the momentum space index. In absence of the magnetic field H the bands for 
the two spin species are degenerate. In the special case of if = and = es, up and 
down spins occupy a single band. For H — and ^ ^b, the band gap for each spin 
species equals (es — e^)- At half-filling the ground state is obtained by completely filling 
up the lower bands which leads to the insulating behaviour of the system. Interesting 
features are expected with the application of a magnetic field in the system. Switching on 
of H removes the spin degeneracy of the bands. Upon increasing H from zero the upper 
and lower bands of up-spin go down in energy scale while the corresponding bands of 
down spin go up. As a result of this there appears a critical value of H, He — {sb — e^), 
where the electrons de-populate the down-spin band and go over to the upper vacant 
band of up-spin electrons. This process makes the bands partially filled and one expects 
a metallic behaviour. This aspect must be reflected in the Drude weight (D) which is a 
measure of the dc conductivity. In the next section we discuss the results of the present 
calculation. 

III. Discussion of Results 

We consider the case with = while the parameters e^, U and T could be varied. 
Fig.l shows the Drude weight {D) as a function of the magnetic field (H) at half-filling 
with = (and, therefore, eb — ^a) and U — ior different system sizes for low and 
high temperatures. Let us first consider the low temperature case {T ^ 0). Decrease in D 
with magnetic field H is the characteristic feature of all the three curves corresponding to 
N = 54, 102, 154, which could be anticipated from the energy band equations obtained in 
the previous section. For 63 — 0, U—0 both the species of spins belong to a single band. 
At half- filling for H—0, D is maximum. But with the increase in H the degeneracy 
of the up and down spins is lifted up. The up-spin band goes down and the down-spin 
band goes up. So, the uppermost down spin electrons can "see" vacant lower energy 
levels in the up-spin band and emigrate to occupy these levels. Consequently both the 
bands will be driven away from half-filling, and therefore, D decreases with H. The step- 
like fall and the oscillation in D are finite size effects. Finite system size accounts for 
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the discrete values of momentum index k or, in other words, leads to the discretization 
of the energy levels. As a result of this, de-excitation (i.e. flipping of the down-spins) 
occurs at specific values of H which is responsible for step-like fall in D. This feature 
disappears with increasing system size as can be observed from the plot. For a system of 
infinite size there will be no change in the spectrum with the variation of because the 
spectrum is continuous, there will be merely a change of variable from k to k' {k'~k + (t)i). 
Whereas, for a finite-size system the spectrum is discrete i.e. k takes on discrete values 
{k = 2Tm/L, n = 0, 1, 2, L — 1, where L is the system size). A change in will shift an 
individual energy level to a new value where there was no level existing before. In other 
words, the spectrum oscillates with ^i, which is again refiected in D. This behaviour also 
vanishes with increasing system size. On the other hand, there is a smooth fall in D as 
a function of H for eb — 0,U — at finite temperature {kT — 0.1). The step-like fall 
and the oscillations, as observed in case of very low temperature, are absent here. This 
is because at finite temperatures the higher states above the Fermi level are occupied. 
Therefore, average occupations of up and down spin band change in a continuous fashion 
with the increase in H. Fig.l suggests that the finite-size effects are almost undetectable 
even for a system size A'^ 150 near zero temperature while such effects are smeared out 
much more rapidly at high temperatures. 

Fig.2 shows the plot of D vs. H for = 1 and U — at low and high temperatures. 
Cases of N — 54, 102, 154 are studied at kT — 0.01. In this case each of the spin species 
has two bands and at half-filling the lower bands are completely filled. Therefore, D 
remains zero from i/ = up to a certain critical value of H, say He- When H reaches the 
critical value H,. « (e^j — e^), the upper vacant band of up-spin gets populated at the cost 
of depopulation of the down-spin band. This is exactly what we observe in Fig.2. The 
rise in D with H at higher temperature {kT — 0.1), however, is rather gradual. Finite 
occupation of higher levels at finite temperature is responsible for smoothening of the 
rise. Also the value of D is higher than that at T « 0. In Fig.3 we observe the effect of 
electronic correlation U. Introduction of U destroys the CDW ordering induced by the 
chemical modulation. As a result of this the charge ordering energy gap is reduced, which 
in turn accounts for the lowering of the critical value of i7 at T fa 0. At finite temperature, 
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for reasons mentioned earlier, the rise in D is again gradual. Due to thermal excitations 
across the CDW gap the conductivity is expected to increase with T upto H < Hc{T « 0) 
while it would decrease with T for H > H^iT « 0) which already appeared to be metallic 
(D ^ 0) at r « 0. 

Fig.4 shows the variation of I? as a function of U/t for different values of at T 
^ for — 54. In this plot, a fixed strength of the chemical modulation is considered 
(cb = 0.5). It is interesting to note that the rate of change of D is very much dependent 
on the value of the field H . It is really governed by the nature of the ground state as 
decided by the competition between e^, U and H. This system favours a CDW ground 
state in absence of correlation and field {U — Q,H — Q). In this phase D is zero [17]. 
For a small value of H (=0.3, well below He) this CDW ordering is still retained. Now 
the increase in U/t breaks up the CDW ordering and as a result of this the Drude weight 
shows a rise upto a certain U/t = {U/t)c where the SDW instability sets in. In the SDW 
phase the Drude weight will fall in a characteristic fashion; the peak in D really marks 
this transition [17]. To verify this conjecture we plot the CDW order parameter (C), the 
SDW order parameter (S) and the magnetization (M) as functions of U/t in Fig. 5(a), for 
H — 0.3 and cb — 0.5. The CDW order parameter shows a gradual fall with increasing 
U/t while the magnetization remains zero upto {U/t)c- Beyond this value of U/t the SDW 
order parameter takes over (Fig.5(a)). Note the "knee" in the plot of C together with 
a small rise in M associated with the CDW/SDW transition. For a value of H ^ He, 
D ioT U/t — because of depopulation of down-spin band as mentioned earlier. 
Hence, the rate of increase of D with f7/t is of much appreciable magnitude (the case of 
H = 0.5 in Fig.4). This is because, at this value of H, the overlap between the lower 
band of the down spins and the upper band of the up spins increases. Therefore, with 
the increase in U/t, the double occupancies are broken and simultaneous depopulation 
of the down-spin band takes place. This enhances the conductivity. This is reflected in 
Fig.4 {H — 0.5) as well as in the corresponding plots of C and M as functions of U/t in 
Fig.5(b). It can be seen from Fig.4 that D saturates to a high value over a small range 
of U/t. This region on the U/t axis essentially corresponds to a magnetically polarized 
phase as revealed in Fig.5(b) which shows comparable values of C and M in this regime. 
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Much higher values of U/t correspond to an SDW phase where the Drude weight falls 
sharply [H — 0.5 in Fig.4) with a fast rise in S (Fig.5(b)). However, for very large values 
of H (> He), the depopulation of the down-spin band is nearly complete. So very small 
changes in C and M (as functions of U/t) are expected in the CDW region and this is 
really observed in Fig. 5(c). Consequently the Drude weight shows a very small increase in 
its value as plotted against U/t in Fig.4 (the case of H — 0.8). Here, in Fig. 5(c), we can 
identify two values of U/t, say (C//t)ci and {U/t)c2 between which M clearly dominates 
over C indicating a phase of strong magnetic polarization. In the corresponding region in 
Fig.4 [H — 0.8) we find D saturating to its maximum. Beyond {U/t)c2, however, we find 
usual drop in D (Fig.4) together with a sharp rise in S in the SDW sector. We note that 
the changes in the order parameters as well as the Drude weight become sharper with 
larger values of the field. 

IV. Conclusion. 

Summarizing, we have studied the conductivity of finite-sized Hubbard ring threaded 
by a magnetic flux in presence of an alternating modulation in the site potentials. Sys- 
tems at half-filling are considered at low and high temperatures. We have considered 
system sizes ^ 100 and found that the finite size effects in Drude weight disappear for 
systems consisting of nearly 150 sites at low temperatures. The competition between the 
correlation, the magnetic field and the chemical modulation tune the conductivity in an 
interesting fashion. In absence of correlation and modulation in site potentials there is 
a slow fall in the Drude weight as a function of the external field. The oscillations and 
step-like fall in the Drude weight quickly disappears with increasing system size. Upon 
introduction of the chemical modulation there arises a critical value of the magnetic field 
across which there appears a sharp rise in the conductivity. This accounts for a transition 
from an insulating phase to a conducting one. However, the introduction of the Hubbard 
correlation tends to destroy the charge ordering effect induced by the chemical modulation 
and, consequently, reduces the value of the critical field. This accounts for an enhance- 
ment of conductivity as well. The critical field, however, is found to be insensitive to the 
system size. The competition between the three factors axe visible from the variation of 
the SDW and CDW order parameters and the magnetization of the system. Transition 
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from a CDW to an SDW phase (at low field) or that from a CDW to a magnetically polar- 
ized one to an SDW (at high field) are found to have notable effects on the conductivity 
of the system. At finite temperatures the nature of variation of the conductivity with the 
external field is much similar to that observed for T but the rise or fall in D is much 
more gradual. Also the finite size effects are not visible even for a system size of A/" 50 
at higher temperatures. It seems from the present study that the interplay of correlation 
and chemical modulation brings in observable features in the magnetoconductance in a 
system. Effects of variation in band filling and the role of diagonal/off-diagonal disorder 
could be interesting areas of further study of this problem. 
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Figure captions: 

Fig.l: Plot of the Drude weight D as a function of magnetic field H for three different 
system sizes: = 54, 102 and 154 at kT = 0.01, and for two different system sizes = 54 
and 102 at A;r=0.1 for = 0.0, U — 0.0. The scale of energy is t — 1.0. is chosen to 
be zero. 

Fig.2: Plot of D versus H for N ^ 54, 102 and 154 at kT = 0.01, and for A^=54 and 102 
at fcT=0.1 for es = 1.0, U — 0.0. The scale of energy is t — 1.0. is chosen to be zero. 
Fig.3: Plot of D versus H for N ^ 54, 102 and 154 at fcT=0.01, and for A^=54 and 102 
at kT=0.1 for = 1-0, U — 1.0. The scale of energy is t = 1.0. is chosen to be zero. 
Fig. 4: Variation of D as a function of ?7/t at a fixed value of chemical modulation e_B = 0.5 
for H = 0.3, 0.5 and 0.8 at kT~0.01. The ring consists of 54 sites, is chosen to be zero. 
Fig. 5: Plots of CDW order parameter C, magnetization M and SDW order parameter S 
against U/t for (a) H = 0.3, (b) H = 0.5 and (c) H = 0.8 at kT = 0.01. The ring consists 
of 54 sites, = 0.5 and = 0. 
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Fig. 2 
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Fig. 3 
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